Mode-Locking and Mode-Competition in a Non-equilibrium Solid-State Condensate 
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A trapped polariton condensate with continuous pumping and decay is analyzed using a general- 
ized Gross-Pitaevskii model. Whereas an equilibrium condensate is characterized by a macroscopic 
occupation of a ground state, here the steady-states take more general forms. Some are characterized 
by a large population in an excited state, and others by large populations in several states. In the 
latter case, the highly-populated states synchronize to a common frequency above a critical density. 
Estimates for the critical density of this synchronization transition are consistent with experiments. 
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Recent experiments [H, [3, [3] have provided substan- 
tial evidence for a new type of Bose-Einstein condensate, 
formed from polaritons in semiconductor microcavities 
Although in many respects these results parallel 
those of condensation in atomic gases, the similarities 
conceal some fundamental differences. In particular, the 
lifetime of a polariton is typically only a few picoseconds, 
and is less than the lifetime of the condensate 0. The 
condensate is therefore a non-equilibrium steady-state, 
in which the decaying polaritons are continually being 
replenished. 

Several consequences of this non-equilibrium aspect of 
the system have now been predicted, based on both mi- 
croscopic calculations and generalizations of the Gross- 
Pitaevskii approach [3, 13 0> HI • While these predictions 
are undoubtedly interesting, more dramatic departures 
from the physics of equilibrium condensates are seen ex- 
perimentally. An equilibrium condensate is characterized 
by a macroscopic occupation at the chemical potential; 
for a trapped, ideal Bose gas this is the ground-state 
energy of the trap, while more generally it is the low- 
est eigenvalue of the Gross-Pitaevskii equation. In some 
experiments there is indeed evidence for equilibrated po- 
lariton distributions and large occupations of the k = 
ground-state of untrapped polaritons [l], [3, [3] , but in oth- 
ers there are massive occupations associated with excited 
states. This has been observed in pillar traps [9(, lattices 
[13], and disorder 11, 12 1. Such results have been ar- 
gued to be evidence of polariton lasing [3, |T3] , distinct 
from condensation, but this distinction is clear only to 
the extent that condensation is restricted to equilibrium. 

The aim of this paper is to outline a theory of the 
trapped polariton condensate with pumping and decay. 
The treatment is based on a recently-proposed general- 
ized Gross-Pitaevskii equation (gGPE) [3,13] ■ Motivated 
by the recent exp eriments, we go beyond the assumption 
[3, 13, 0, 13, 0, El of condensation at a single ground-state 
energy, and systematically investigate how the dynamics 
of the pumped system can lead to more general forms. 
We first consider a small trap at low densities, and show 
that the gGPE reduces to kinetic equations for the occu- 
pations of the trap states. The occupations of the states 



are determined by gain and loss processes which are in- 
dependent of energy, i.e. a laser-like mode competition. 
Thus the steady-state solutions can include massive oc- 
cupations of an excited state, or of several states simul- 
taneously; the character of the steady-states in this limit 
can be predicted from the wavefunctions of the trap. 

Having established the steady-state structure at low 
densities, we then numerically explore how it changes 
away from this limit. We demonstrate that, if the trap is 
such that there are several massively-occupied states at 
low densities, then above a threshold density the steady- 
state reverts to a massive occupation at a single energy. 
This can be understood as the classical phenomenon of 
synchronization (mode-locking) of coupled nonlinear os- 
cillators [161 ] . Supposing that the dominant nonlinearity 
is the polariton-polariton interaction we predict that, for 
realistic parameters, polariton systems could be tuned 
through the synchronization transition. They thus give 
access to both an interaction-dominated regime, and a 
regime of laser-like mode competition. 

Synchronization in the polariton condensate has been 
independently considered by Wouters [l7[. That treat- 
ment is complementary to this one, considering the spe- 
cific problem of a double- well trap using a related model. 
The double- well has also been considered under resonant 
pumping [13] , which further differs from the present work 
because the condensate is directly induced by the pump. 

We consider the order-parameter equation [3] 

id t il> = [H + U\i>\ 2 + i(7eff(r) - r|V| 2 )] V- (1) 

This complex Ginzburg-Landau form [13] has been ar- 
gued [3, [3] to be the generic mean-field description of a 
polariton condensate with continuous pumping and de- 
cay. The first two terms on the right form the standard 
Gross-Pitaevskii mean-field equation, with local density 
\ip\ 2 . Hq includes a quadratic approximation to the po- 
lariton dispersion, and a potential due to disorder or the 
trap. U > is the strength of the repulsive interaction 
between polaritons, which is treated as a contact interac- 
tion because its range is of the order of the exciton Bohr 
radius. The polarization degeneracy of the polaritons has 
been neglected for simplicity; it could be treated using a 
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two-component order-parameter equation [141 ]. 

The last two terms on the right of (QJ account for the 
pumping and decay. The pumping model involves a reser- 
voir of high-energy particles, created by some external 
excitation. The condensate is populated by stimulated 
scattering from this reservoir, contributing a linear gain 
term i~/(r)?p, where 7(7-) is related to the reservoir den- 
sity. This term combines with a similar term from the de- 
cay of the polaritons, giving the overall linear gain term, 
with coefficient 7 c ff , in |T]). However, the rate of conden- 
sate growth should reduce with increasing density, as the 
pump reservoir becomes depleted. This effect, modeled 
using a single gain-saturation coefficient T, gives the fi- 
nal term in ([1]). Physically, this form of pumping can 
be interpreted in terms of a pump which tries to locally 
enforce a steady-state density 7 e g(r)/r. 

Strong- Trapping Limit- The steady-states of ([1]) can 
be determined analytically in the limit of strong trap- 
ping, where the nonlincarities are weak compared with 
the single-particle level spacing. We may then treat them 
with degenerate perturbation theory. We expand t) 
in terms of the eigenstates of H$, ip^, and retain only 
resonant terms in the resulting equations. To simplify 
the notation we analyze a trap with only two single- 
particle states, so tp(r,t) = ^\{t)ipi(r) + /^(i)^^), anci 
assume homogeneous pumping. The amplitude Hiit) 
{pv{t)) obeys 

*Ai(2) = [ E i(2) + *7cff + (U- iT) 

x (^1(2) |^i(2) I 2 + 2&| M2(1) | 2 )] ^ 1(2) . (2) 

Ei is the single-particle energy, and the wavefunctions 
have been taken to be normalized and real. 771 = 
""1111,772 = 1*2222, arid 6 = 1*1122 are matrix elements for 
a local nonlinearity, 



(3) 



They parametrize the inhomogeneous density profile of 
the trap states, with 771 and 772 describing the inhomo- 
geneity of the states, and 6 their overlap. 

Introducing number and phase variables in a rotating 
frame, 



Mi(2) (t) = e-W* ^/n 1(2) (t)e-^w {t) , (4) 

separates the number and phase dynamics. The former 
obeys the rate equations 

7J 1(2 ) = 2rni (2 )(7cff/r - 771(2^1(2) - 26n 2 (i)), (5) 

with terms describing the stimulated scattering from the 
reservoir and the spontaneous decay, and the reservoir 
depletion. This result can be understood as a gener- 
alization of the kinetic description of polariton lasing 
to treat the spatial structure of the trap. 



It describes the extensive component of the occupation, 
and hence effects such as spontaneous pumping are miss- 
ing. They would be important for a finite system close 
to threshold [23 |. 

The phase dynamics is straightforward, obeying 

0i(2) = E 1{2) - C/(7cff/r - ?7i(2)ni(2) ~ 26n 2 (i)). (6) 

Each mode oscillates freely, at a single-particle energy 
which is shifted by the repulsive interactions. 

The rate equations ([5]) have several steady-state solu- 
tions above the bulk condensation threshold, j e g = 0. 
There are always two solutions corresponding to conden- 
sation in each of the trap states: either n\ = ^^/(riiY) is 
finite and tt^ = vanishes (state Si), or vice versa (state 
S2). However, if either 771, 772 > 26 or 771, 772 < 2b there is 
also a steady-state T, with massive occupations of both 
trap states: 

26- 

n±(2) 
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The left panel of Fig. [T] shows the parameter regions 
in which the different steady-states are stable. Lineariz- 
ing ([5]) about the steady-state Si, we find that popu- 
lation fluctuations decay with rates Ai = 27 and A2 = 
27(26/771 — 1). Thus this solution is stable for 771 < 26. 
This is the condition that the occupation of the first trap 
state alone, determined by the pumping and the self-gain- 
saturation parameter 771 , is sufficient to keep the second 
below threshold. The analogous condition 772 < 26 holds 
for the existence of a stable solution in which only the 
second trap state is occupied, S 2 - If neither criterion is 
satisfied both states must be occupied, and this is where 
the two-mode solution T is stable - fluctuations there de- 
cay with rates 27 and 27(26 — 771) (26 — ?72)/(46 2 — 771772). 

The different steady-state solutions can be realized 
in many different potentials. For two widely-separated 
states we generically have 6 — > 0, and hence obtain two 
coexisting condensates. While this result is perhaps ex- 
pected for separated traps, it can nonetheless also occur 
when there is substantial spatial overlap. A simple model 
which demonstrates this is a one-dimensional hard-wall 
trap of size a, with a finite trap of size 6 and depth V at 
its center. The right panel of Fig. Q] shows the different 
steady-state regions for the two lowest eigenstates of this 
potential. 

Experimentally, the presence of massively-occupied 
states in polariton systems is shown by bright lumines- 
cence peaks. The characteristic frequencies in such op- 
tical spectra follow from ((6|). In Si(S 2 ) the massively- 
occupied states will lead to strong emission at the fre- 
quency Ujcs/T + Ei(n), while in T both peaks appear 
simultaneously. Note that here the blueshift of the con- 
densing modes is Uj c s/T, irrespective of which steady- 
state we consider. Physically, this is because the energy 
shifts are determined by the density, which is fixed by 
the pumping. 



FIG. 1: Left panel: steady-state structure for condensation in 
a two-state trap, for the generalized Gross-Pitaevskii model 
neglecting non-resonant terms ([2]| . Lettering denotes the sta- 
ble steady-states in each region, with condensation in both 
trap states (T), state 1 alone (Si), or state 2 alone (S2). In the 
lower-left region both Si and S2 are stable, and the steady- 
state is selected by the initial conditions. Right panel: steady- 
state structure for the two lowest states of a one- dimensional 
model trap, consisting of a hard-wall trap of length a with a 
trap of length b and depth V at its center. 



In addition to the strong emission associated with the 
condensing states, we also expect peaks in the optical re- 
sponse associated with the non-condensing trap states. 
For the mean-field model in the resonant approximation 
these states are not populated, and hence would appear 
in the absorption but not the luminescence. However, 
they could develop non-macroscopic populations due to 
effects beyond that model, in which case they would ap- 
pear weakly in luminescence. These peaks are shifted by 
the mean-field repulsion with the condensate, so in Si 
the absorption peak will be at E2 + 2bn\. Since its width 
will be proportional to A2, the mode is narrow close to 
the critical line r\\ = 2b, where it goes unstable. 

Beyond Strong- Trapping- Having established the 
physics of ([T]) in the strong-trapping regime, we now in- 
vestigate how it develops with increasing nonlinearity. 
To do this we use real-space discretization to solve ([1]) 
directly, for the one-dimensional model described above. 
This potential is chosen as a simple, few-parameter 
model, which can be tuned into the different steady-state 
regions. We choose a = m = T = 1, taking a as the unit 
of length, H 2 /(ma 2 ) as the unit of energy, and working 
with a rescaled density \ip\ 2 — > \ip\ 2 /F. Both the repul- 
sive interaction U and the gain saturation T are believed 
to be significant for polaritons [3, 0, @] , and so for this 
demonstration we take U/T = 1. 

Fig. [2] illustrates how the steady-state spectra develop 
with increasing pumping. The analysis above, applied 
to the two lowest states of this potential, predicts the 
two-mode steady-state T. This is in agreement with 
the results at weak pumping (lowest panel). Increas- 
ing the pumping there is the overall blueshift associated 
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FIG. 2: Spectral analysis of the polariton field in the steady- 
states of the generalized Gross-Pitaevskii model, with the po- 
tential described in the caption of Fig. [1] The pump strength 
increases through 7 e g = 1, 30, 60 from the lowest panel to the 
highest. The grayscale is the computed amplitude \ij)(E,x)\, 
normalized by x/TWr to account for the overall increase in den- 
sity. U/F = l,b = 0.05 and V = 25. 

with the increased density. The lowest two modes still 
dominate the spectrum, but their energy splitting has 
reduced slightly, and several further emission peaks ap- 
pear (middle panel). Further increasing the pumping, 
the spectrum switches to emission at a single frequency 
(top panel). 

At a general level, these results are expected conse- 
quences of the non-resonant terms dropped from @. In 
a two-mode model with states of opposite parity, for ex- 
ample, there is the additional Josephson term 

&(E/-ir)/^ (2)J u 2 (i)M2(i) (8) 

in the equation j2|) for ;Ui(/X2). Thus the condensing 
states can drive nonlinear emission at other frequencies, 
as we see at intermediate pumping. Furthermore, the 
non-resonant terms couple together the phases of the con- 
densing modes in ((6]) . Thus we expect frequency pulling, 
and eventually synchronization of the condensing states 

Coupled oscillators synchronize to a common fre- 
quency when the phase-phase couplings become compa- 
rable to their energy splitting A. Since here the states 
have similar and overlapping density profiles, the matrix 
elements ([3]) are all of order one. The scale of the phase- 
dependent couplings in ^ is then set only by the non- 
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linearity \JU 2 + T 2 and the polariton density p « j e s/T. 
We may thus estimate the critical polariton density for 
the mode- locking transition. p, : . from 

A ~ Pc ^u 2 + r 2 ~ 7crit v/i + c/ 2 /r 2 . (9) 

This form is consistent with Fig. [21 and with simulation 
results (not shown) for other values of the nonlinearities 
U and r. In particular, in the same model the synchro- 
nization occurs between 7 c ff/r = 20 and 40 when [7 = 0, 
and between U/T — 20 and 40 at fixed 7 e ff/r = 1. 

Synchronization due to nonlinear gain is well-known 
as the basis for mode-locked lasers. However, an impor- 
tant difference between the polariton condensate and a 
laser is the presence of strong interactions U between the 
particles, due to the excitonic component of the polari- 
ton. Thus frequency-pulling and synchronization could 
be expected to have a much wider role in the physics of 
polariton condensates than they do in lasers, occurring 
on large energy scales at low intensities. 

Let us estimate the critical polariton density p c for 
synchronization, supposing it is controlled by the real 
nonlinearity U. Estimates of this interaction are available 
for plane-wave excitons in a perfect quantum-well, and 
the localized exciton states of a disordered quantum well 
[25l | . For current experiments, a plausible upper limit is 
the result U ~ Q,n{m x W)~ 1 from the disordered models, 
with disorder energy-scale W ~ lmeV. m x ~ 0.5m e is 
the exciton mass, and f2# ~ 20meV the Rabi splitting. 
Thus the phase boundary for synchronization in a trap 
of scale L t is 

n R {m x W)- 1 Pc ~ l/(mi 2 ). (10) 

For a trap of L t ~ 5/im [9j and a polariton mass m ~ 
10" 5 m e this gives p c ~ 10 n cm~ 2 , or n c ~ L 2 p c ~ 10 3 . 
This estimate is about one order of magnitude larger than 
the densities usually reported for polariton condensates. 
Thus, with a suitable potential, coexisting polariton con- 
densates of different frequencies should be expected in 
tight traps. The synchronization transition could then 
be observed with increasing density, if necessary using 
larger traps to reduce p c into the experimentally accessi- 
ble range. 

Since we do not expect synchronization at current 
densities in tight traps, the present theory is consistent 
with the observations there of well-resolved emission lines 
HH • It also appears to be consistent with the existence 
of long-range coherence [3, EH at current densities, with 
softer traps formed by disorder. 

Concluding Remarks- Because the tightly-trapped po- 
lariton condensate has well-resolved emission lines, it 
could provide a sensitive probe of the physics of non- 
equilibrium condensates. In particular, the present re- 
sults will allow the order-parameter equation |TJ) to be 
tested. While this form certainly captures much of the 
physics [B|, Q , effects which are important elsewhere are 



missing. Most obviously, there is no thermalization with 
the reservoirs [26|, [23] , which exchange particles with the 
system irrespective of energy, and do not directly cause 
transitions between trap states. Such effects would ap- 
pear in generalizations of the kinetic equations ([SJ • Since 
some groups [HQ report thermalized distributions, such 
generalizations may prove necessary. 

To conclude, we have considered the trapped polari- 
ton condensate in the framework of a generalized Gross- 
Pitaevskii model. At low densities, this model admits 
solutions which, differently from an equilibrium conden- 
sate, involve massive occupations of excited states, or 
of several states simultaneously. We have derived cri- 
teria for predicting the nature of the steady-states in 
a given geometry, and shown that the steady-states are 
selected by gain-competition effects. At a general level 
such physics is of course familiar in lasers, though it has 
not previously been considered for the polariton system. 
Moreover, the direct interactions between polaritons cre- 
ate differences compared with the photon laser: blucshift- 
ing the modes in the weak-nonlinearity regime, and caus- 
ing frequency-pulling and synchronization at stronger 
nonlinearities. Our estimates for the critical density of 
the synchronization transition suggest that it could be 
cleanly observed in tight traps [9(, and may be respon- 
sible for the observations of long-range coherence across 
disorder potentials 0, [ll| . 
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